Asymptotic symmetries in an optical lattice 
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It was recently remarked by Lutz [Phys. Rev. A 67 (2003), 051402(R)] that the equation for the 
marginal Wigner distribution in an optical lattice admits a scale-free distribution corresponding to 
Tsallis statistics. Here we show that this distribution is invariant under an asymptotic symmetry of 
the equation, hence that this scale-free behavior can be understood in terms of symmetry analysis. 
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In a recent paper with R. Mancinelli [1] we have consid- 
ered asymptotic symmetries of differential equations and 
of their (asymptotic) solutions; in there we considered 
in particular appUcation to a given class of anomalous 
reaction-diffusion equations which had been extensively 
studied numerically [2], obtaining a theoretical explana- 
tion of the observed long-time behavior of solutions. The 
same approach has also been extended to a discrete ver- 
sion of these equations [3] . 

In this note we apply our approach to a different kind 
of anomalous diffusion; that is, we focus on the equation 
describing anomalous transport in an optical lattice [4, 5] . 
The equation for the marginal Wigner distribution w(p, t) 
of the momentum p at time t reads 



d 

wt = [h{p)w - g{p)wp] 

where the functions h{p) and g{p) are given by 



Hp) := 
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1 + {p/pcY 



9{p) ■■= 70 + 
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1 + {pIpc? 



(1) 
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here a, 70, 71 are certain constants, Pc represents the cap- 
ture momentum. 

Equation (1) should be complemented by a choice of 
the function space to which the function w{p,t) belongs. 
On physical grounds, we require that w{p, t) at fixed t is 
normalizable and vanishes asymptotically. 



\'w{p, t)| ^ for \p\ 



00 
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We stress that there is a range of parameters (—1/2 < /i) 
[5] for which there is an equilibrium distribution with 
these properties; this is the range to be studied here. 

The explicit expressions of h(p) and g{p) make that, 
with the shorthand notation P{p) := 1/[1 + (p/pc)'^], the 
equation (1), (2) reads 



Wt 



(70 + /3(p)7i) Wpp + 
+ {a(3{p) - 27i[/3(p)/p,]2) pw. 
+ a {(3{p) - 2[{p/p,)P{p)]^ 



(4) 
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Attention to eq. (4) was recently called by Lutz [5] 
(to which the reader is also referred for derivation and a 



discussion of this equation), who remarked that - quite 
surprisingly - the equilibrium distribution is a Tsallis 
i.e. a distribution optimizing the Tsallis 
Indeed, the stationary solution of (4) turns 



distribution, 
entropy [6]. 
out to be 
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where 



f3 



a 



2(70+71) 



270 
apl 



(5) 



1 + S, (6) 



and Z is a normalization factor, which can be chosen so 
that J_^wo{p)dp = 1. The physical range is q < 3; 
the case 5/3 < g < 3 (for which the second moment 
/ p^wo{p)dp diverges) corresponds to anomalous moment 
diffusion [5]. 

The peculiar properties of this equation and its sta- 
tionary solution, in particular concerning the power-law 
decay of wq{p) for large \p\, hence its scale-free nature, 
has been studied by Abe [7] in connection with dilation 
symmetries and canonical formalism. 

Here we will consider generalized scaling transforma- 
tions; that is, local transformations acting as a standard 
scaling in the independent p and t variables, and as a p 
and t dependent scaling in the dependent variable w{p, t). 
Moreover, we will not require to have a full invariance of 
the equation, but be satisfied in asymptotic invariance 
(for large \p\ and t and hence, in view of (3), small w), 
as discussed in detail in [1]. 

We start by recalling a general feature of symmetry 
transformations and symmetry invariance for PDEs and 
their solutions [8-10]; we will specialize to a scalar PDEs 
for w = w{p, t). 

Under an infinitesimal transformation with generator 



X 



^{p, t, iv)dp + t{p, t, w)dt + (flip, t, w)du 



(7) 



(that is, p p = p + s^{p, t, w), and so on) the function 
'w{p, t) is mapped to a new function w(jp, t) with w{p, t) — 
'w{p, t) + £[6'w{p, t)] and 



Sw = [(p — Wp^ 



WtT\, 



v{p.t) 



(8) 



Applying this to the distribution (5) and to generalized 
scalings we get with easy computations that (up to a 
common factor) the transformation generated by (7) with 
^ = —p leaves (5) invariant if and only if 
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obviously t = at remains unrestricted at this stage. 
Conversely, a function w{p) is invariant under X = 
-pdp + {2f3Zi-'^)p'^w''d.^ if and only if w{p) is of the 
form (5); this is easily checked by using (8). We will 
write, from now on. 



1 , z/ = 2I3Z^ 



(10) 



Albeit our analysis will be conducted at the infinitesi- 
mal level, we note that the action of the one-parameter 

group (denote by s the group parameter) generated by 
X = —pdp + yp^w^^^d^ on the p and w variables is to 
map {p,w) = {po,wo) into (p(s),w(s)) with 

p{s) = e-^po , 

w{s) = [1 + {upl/2){l - e-'^w'o] wo. ^ ^ 

We are now going to discuss what is special in the 
generalized scaling vector fields identified above 



X = -pdp + atdt + {up'^wi)du 



(12) 



for what concerns their action on the equation (4). 

We follow the general procedure for group analysis of 
differential equations [8-10]. By acting with X on the 
(p, t, w) variables, we induce a transformation in the par- 
tial derivatives of w with respect to p and t as well (also 
called the second prolongation of X); this is described - 
restricting to partial derivatives occurring in (4) - by the 
vector field 



X + + "^p^ 



d_ 

' dwt 



d 



(13) 



pp 



the coefficients ^ appearing in here is as follows: un- 
der the infinitesimal transformation described by X , the 
partial derivatives get transformed as Wp ^ Wp + s^p, 



Wt 



Wt + e*t, Wp 



When X is the generalized scaling (12), i.e. for the 
mapping p — > (1 — e)p, t ^ (1 + ea)t, w — > (1 + 
e{vp'^w^))w, the corresponding ^ can be easily computed 
by the so-called prolongation formula [8-10] to be 



*t = 

*P = 



—crwt + vqp^w^ 
Wp + 2vpw^^+^'^ + i/(l -I- 5)p'^w^ 
2Wpp + 2vw^'^+^^ + 



(14) 



vw'^^ ^) (1 -I- ^) [Sp^Wp + pw{Awp + pwpp) 



These determine how the differential equation (4) of in- 
terest here is transformed under X. 

In computational terms, this is done by applying Y on 
(4), and substituting for Wt according to (4) itself; the 
obtained expression must be zero for X to be a symmetry 
of the equation. Proceeding in this way, we obtain a 
condition which we write in compact form (see below for 
the explicit expressions of the functions A/-) as 

Ao{p,t,w) + Ai{p,t,w)wp + A2{p,t,w)wpp = 0. (15) 

The functions Ak must vanish separately for (15) to be 
satisfied. 



The explicit expression of A2 is 
A 7ipg(pg(2 + a) + (4 + a)p^) 

A2 = -7o(2-hcr) r , 010 • (lo) 



[Pc +P^]^ 

The limit of this for large \p\ is nonzero unless we choose, 
as we do from now on, 



a = -2 . 
With this choice, X reads 

X = -pdp - 2tdt + up^w'^du 
and A2 reduces to 

-271 Pc^p2 



Ao = 



(Pc^+p'^) 



2 ' 



(17) 



(18) 



(19) 



therefore, A2 ^ for \p\ 00. 

As for Ai, with the choice (17) it reads 

^1 = [2p{{pl+p''){api-2{l + S)'yoiyw^{pl+w^)^) 

-27ip2((l + S)vp'^.w^ +P^(-1 + vw^p'^ + Svw^p^) 
+pI{1 + 2vw^p' + 25y'w^p')))] x [pl+p"]'^ . 

(20) 

This is a more involved expression, but it is still easy to 
check that in the limit \p\ — » 00 and w — *■ 0, recall (3), 
we get A\ 0. 

Finally, and again with (17), Aq reads 

^0 = [w{-2yw^{pl+p'^){-fipl{pl-p^)+ 

+7o(Pc+p')')+ 

+apl{-{{2 + 6)vw''p^) - 2p2p2(3 _^ 2vw^p'^)+ 
+pI{2 - 2vw^p^ + (51/wV)))] X 

(21) 

Once again, for \p\ — > 00 and «; ^ we get ^0 — * 0. 

This computation shows on the one hand that the one- 
parameter group of transformations generated by (18) is 
not a symmetry of the equation (4), as the A^ are not 
identically zero; on the other hand, it shows that (18) 
is an asymptotic symmetry for this equation, defined in 
the fimction space identified by (3). No vector field of 
the form (12) with u ^ — 2 is an asymptotic symmetry 
of (4); hence (18) is the only invariance generator for (5) 
which is also an asymptotic symmetry of (4). 

Let us discuss the effect on (18) on the normalization 
condition (3). Using (8) we obtain that a generic function 
w{p,t) is mapped by (18) into w = w + eSw with Sw = 
vp^w'^ +pwp + 2twt] on solutions to (1) this reads 

d 

Sw = up^w'' + pwp — 2t —{hw — gwp) . (22) 

With I[w] := w{p,t)dp, we obviously have I[w] = 
I[w] + el[6w]. In view of (22), the latter amounts to 

I[Sw] = ullp^w'^] + I{pwp] +2tl[dp{hw-gwp)] . (23) 
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It is immediate to sec, using also an integration by parts 
for the last one, that (3) implies the finitcness of the last 
two integrals in (23). 

As for /[p^w^], note that (3) holds if, for \p\ oo, 
w{p) ~ l/p^*^ with k > 1/2. For general solutions w{p, t), 
the condition of normalizable variation kq > 3/2 is more 
restrictive than the normalization condition k > 1/2 for 
all q in the physical range q < 3. 

Note that for wip) given by (5), k = 1/S and the condi- 
tion kq > 3/2 holds for all g < 3; that is, functions which 
are near to the stationary solution wo(p) will always have 
normalizable variation under (18). 

Thus, strictly speaking, our method should be applied 
only on functions satisfying |/[p^w'^]| < cxo, i.e. decaying 
(for |p| oo) faster than l/|p|^/'. As observed above, 
this includes all functions near to the Tsallis distribution 
(5), for all q in the physical range. 

The (asymptotic) invariance of (4) under (18) could 
also be analyzed using the systematic procedure of [1]; 
this requires to introduce symmetry-adapted coordinates 
{v,y,a), with v the dependent variable. These are the 
coordinates in which (18) is simply X = —2cFda, and 
its action on derivatives up to order two is described by 
Y = X + 2v„{d / dvcr) [8-10]. In the present case, adapted 
coordinates are a = t, y = p^ /t, v = — {v5 /i)p'^ . 

Finally, let us briefly mention differences with the 
work by Abe [7]. Abe considered a dilation symme- 
try of (4), based on an auxiliary field K{p) which could 
be not normalizable. The Abe symmetry has generator 
G = J p^p w dp. In the language of symmetry theory [8- 
10] this is a nonlocal symmetry (it depends on an integral 
over p rather than just on the value of w at the point p), 
whereas here we considered local ones. Moreover Abe 
worked in canonical formalism, which requires to select 
a symplectic structure, whereas here we did not consider 



any additional structure. The two approaches are thus 
quite different; it appears that symmetries considered by 
Abe are more general, while the one considered here has 
the advantage of being a local one. 

Let us summarize our discussion. We have considered 
the equation for the marginal Wigner distribution tii{p, t) 
of the momentum p at time t in an optical lattice from the 
point of view of symmetry analysis. It is known that this 
equation admits the distribution Wo{p) as a stationary 
solution. 

We have identified the generalized scaling invariance 
group of wo{p), described by (11) and generated hy X = 
-pdp + {2l3Z'i-^p'w'i)d^. 

We passed then to consider generalized scalings (12) 

acting on the time variable as well; it is easily seen that 
these cannot be symmetries of the equation (4). We in- 
vestigated then if they may be asymptotic symmetries for 
that equation, and observed that a necessary condition 
for this is provided by (17). With this choice, it turns 
out that we have indeed an asymptotic symmetry, whose 
generator is (18), of our equation. 

This shows that the Tsallis distribution (4) is an in- 
variant function for a transformation X which is an 
asymptotic symmetry (for large |p|) of the equation under 
study, and hence that the scale-free asymptotic behavior 

of wo{p) is a, consequence of the asymptotic symmetry 
properties of the equation (4) describing the characteris- 
tics of the optical lattice. 
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